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We discuss the unitary quantum dynamics of the Dicke model (spin and oscillator coupled). A 
suitable quasiprobabilty representing the quantum state turns out to obey a Fokker-Planck equa- 
tion, with drift terms representing the underlying classical Hamiltonian flow and diffusion terms 
describing quantum fluctuations. We show (by projecting the dynamics onto a co-moving Poincare 
section) how the interplay of deterministic drift and quantum diffusion generates equilibration to the 
microcanonical density, under conditions of global classical chaos. The pertinent photon statistics 
reveals macroscopic quantum fluctuations. 



I. INTRODUCTION 

The Dicke model [T] has long been known for its quan- 
tum phase transition (to a superradiant phase as the cou- 
pling exceeds a critical value) [5] and classical as well 
as quantum chaos New interest in the model has 

recently been stimulated by an experimental realization 
(employing a double Bose-Einstein condensate coupled to 
an electromagnetic mode of a surrounding cavity) where 
the zero-temperature phase transition was observed [6j. 
Fascinatingly, a many-body system here exhibits macro- 
scopic quantum behavior involving only two coupled de- 
grees of freedom, one each for spin and oscillator. 

The observation of chaos requires a combination of not 
too small coupling and sufficient excitation. The experi- 
mental realization of such regimes would reveal two dis- 
tinctive signatures of chaos, dissipationless equilibration 
on the energy shell [7] and macroscopic stationary fluc- 
tuations of the photon number, and both of these should 
be accessible in the condensate/cavity setting of Ref. [BJ. 
For large values of the quantum number j fixing the con- 
served square of the spin as J 2 = j(j + 1) >• 1, the quan- 
tum dynamics will be strongly influenced by the classi- 
cal Hamiltonian flow. The essential physics of the Dicke 
model will thus be governed by an interplay of classical 
chaos and semiclassical quantum fluctuations. 

The purpose of this paper is threefold: 

i) We will present a comprehensive discussion of chaos 
in the system. Unlike with the pioneering study 
[5] , we will not limit ourselves to the vicinity of the 
system's ground states but address the full phase 
space available to chaotic fluctuations. 

ii) We discuss how a conspiracy of quantum fluctu- 
ations and nonlinear classical dynamics generates 
'thermalization' to a microcanonical distribution. 
We thus present a first case study where the gen- 
eral role of quantum fluctuations in a chaotic ther- 
malization process is resolved in microscopic detail. 

iii) We will analyse giant fluctuations of observables 
such as the photon number resulting from the uni- 
form coverage of the system's energy shell. 



i) For the classical motion of the Dicke model, we will 
identify regimes of integrable, mixed and fully chaotic 
dynamics in regimes of small, intermediate and strong 
spin-oscillator coupling. Specifically, we will show that 
chaotic trajectories sweep out large intervals of the clas- 
sical action variable corresponding to the number of os- 
cillator quanta. The tendency to uniform coverage of 
the energy shell, typical for dominantly chaotic behav- 
ior, arises already for rather moderate energies, provided 
the coupling is sufficiently large. We conclude the clas- 
sical discussion by demonstrating that the classical Li- 
ouville equation would effectively entail equilibration to 
the energy shell. The disclaimer 'effectively' refers to the 
fact that the classical Q evolves into an infinitely rugged 
distribution of alternating high and zero phase space den- 
sity. Only in the infinite time limit, the ensuing structure 
becomes infinitely filigree, and looks homogeneous at any 
finite resolution. 

ii) To decribe the quantum dynamics, we represent the 
density operator by a suitable quasi-probability, the 
Glauber Q-function aka Husimi function. Schrddinger's 
equation then turns into a Fokker-Planck equation for Q. 
The drift terms therein (first-order derivatives) reflect the 
classical Hamiltonian flow while the second-order deriva- 
tive terms, small in the parameter 1/j, describe quan- 
tum diffusion. The important role of quantum diffusion 
is, under conditions of global classical chaos, to limit 
the shrinking of phase space structures along the clas- 
sically stable directions to a scale of the order j~ 1//2 . 
The competition of classical vs quantum contraction and 
expansion of phase space structures corroborates equili- 
bration to the microcanonical distribution. In particular 
it implies that a) any quasi-probability Q will eventu- 
ally cover the compact energy shell of the model, on a 
time scale bounded from above by the Ehrenfest time. 
Further, our theory reveals explicitly how quantum dif- 
fusion smoothens the ruggedness of classical structures 
over scales ~ j^ 1 ^ 2 - As a result, we obtain thermaliza- 
tion into a genuinely uniform distribution (to be distin- 
guished from a fissured structure that only looks like a 
distribution upon decreasing the level of resolution) on 
time scales of the order of the Ehrenfest time ~ ln(j). 
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It is worth emphasizing that the smoothing mechanism 
relies on the effect of quantum fluctuations during the dy- 
namical evolution, rather than just reflecting the initial 
quantum uncertainty intrinsic to any starting distribu- 
tion. 

iii) The ensuing uniform coverage of the energy shell en- 
tails large fluctuations of observables such as the photon 
number or spin projections. The underlying (semiclassi- 
cal) physics is that during their dynamical evolution, sys- 
tem trajectories fully explore the available phase volume 
in an essentially unpredictable manner. For instance, 
in the superradiant regime and at energies E > huioj, 
where ujq is the spin precession frequency, the average 
photon number is roughly given by n ~ 0(E/huj). The 
fluctuations, 5n superimposed on this mean are of order 
6n~ Oij 1 / 2 ^) ~ 0{n), i.e. 'macroscopic' fluctuations 
as big as the average may occur. 

The rest of the paper is organized as follows. In Sect. [II] 
we discuss the Hamiltonian and the quantum dynamics 
it generates. Specifically, the quasiprobability Q will be 
shown to obey a Fokker-Planck equation. Sect . |III| is 
devoted to the classical dynamics and in Sect. we 
investigate quantum diffusion. Finally, in Sect.^V we 
discuss fluctuations of the photon number by calculat- 
ing the microcanonical averages of moments. In Sect. |VI| 
we point out why the chaotic drift/diffusion mechanism 
bringing about equilibration and large stationary fluctua- 
tions works in other systems as well. Prominent examples 
are the kicked top (recently realized experimentally) and 
the Bose-Hubbard model (a genuine many-body system 
of much current interest). Appendices will detail some 
calculations. 



II. MODEL AND QUANTUM EVOLUTION 

The Hamiltonian of the Dicke model can be written as 
H = h^j J z +uja}a + g^(a + a))J x ^ . (1) 

Here, the operators J a ,a = x,y,z act in a spin-j repre- 
sentation and obey the standard commutation relations 
[J a ,Jb] = \e abc J c , where {e abc } is the fully antisymmet- 
ric tensor. The photon annihilation and creation opera- 
tors fulfill the Bose commutation rules [a, a*] = 1. The 
first two terms in ([!]) respectively describe spin precession 
about the J z -axis with frequency ujq and harmonic oscil- 
lation with frequency u. The last term accounts for spin 
precession about the J^-axis with a 'frequency' oc a + a) 
and for driving of the oscillator by a 'force' oc J x . The 
coupling constant g is a (Rabi) frequency independent of 
H. The appearance of the spin quantum number j in the 
interaction part is owed to the use of the operators a, a' , J 
which arc rather non-classical in character [8]. We note 
that the Hamiltonian ([I]) contains the antiresonant terms 
J + a) + J_a. Only the parity P = expi7r(a^a + J z ) thus 
remains as a symmetry. If the antiresonant terms were 



dropped ('rotating wave approximation'), conservation of 
a) a + J z and thus intcgrability would result. 



A. Coherent state representation 

We aim to explore the quantum dynamics generated by 
the Hamiltonian ([I]) . In view of the largeness of the spin, 
j > 1, we find it convenient to employ coherent states 
which are optimally suited to taking semiclassical limits. 
Specifically, spin coherent states [9HTT] are defined as 
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where \j,j) is a 'maximum- weight' eigenstate of J z , i.e. 
Jz\j,j) = j\j>j)' The states \z) yield the mean values 
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are the three components of a unit vector 
1 = (sin# cos (f>, sin0 sin 4>, cos 0) T , whose angular orien- 
tation is defined through z — e 1 ^ tan(0/2). Higher mo- 
ments reveal minimum angular uncertainty, character- 
ized by the solid angle 4ir/(2j + 1) which defines a 
Planck cell on the unit sphere. The set {\z)} is over- 
complete and allows the resolution of unity by 1 = 
■(l+tz*) 2 \ z ) i z \' T w0 s P m c °h eren t states have the 
overlap (z\z') — ( - 1+z - ^ 



(l+|*PK(l+|z'| 2 )i- 

Similarly, oscillator coherent states are defined as 
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with a € C a complex amplitude and |0) the vacuum, 
a|0) = 0. The resolution of unity in terms of the 
(over)complete set {\a}} reads t=±J da\a)(a\. The 
state | a) assigns a minimal uncertainty product to dis- 
placement and momentum such that these quantities are 
'confined' to a single Planck cell. The latter property is 
also evident from the overlap |(a|a')| 2 = e~\ a ~ a I . 

We shall discuss the dynamics of the system in terms 
of its density operator p(t) and represent the latter by 
the Glauber Q or Husimi function |llj 



Q(a,z) 
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Obviously, that function is real and non-negative every- 
where and exists for any density operator p. By invoking 
the completeness relations given above, one immeadiately 
checks that expectation values of (anti-normal ordered) 
operators of the oscillator can be computed as 
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Expectation values of spin operators can be computed 
by analogous averaging over z- valued functions. We just 
note the example 

(a"V™J a ) = J dadza m a* n (j + l)l a (z,z*)Q(a 1 z) , (7) 

where the unit increment over j in the factor (J + 1) is a 
quantum correction. At any rate, the foregoing proper- 
ties of Q allow us to speak of a quasi-probability density 
which is expected to converge to a classical phase-space 
density in the limit h — > 0. 

To illustrate the use of our quasiprobability we note 
that a density operator projecting onto a pure coherent 
state |zo,c*o) = \ z o) ® | cko) implies a Q-function smeared 



out over a single Planck cell as 

Q(z,a) = ^7^7q~p\2 |(o^K,zo)| 2 ■ (8) 

We shall in fact mostly imagine the system initially pre- 
pared in such a nearly classical state. 

B. Time evolution 

The von Neumann equation dtp = jr[H,p] can be 
rewritten as a partial differential equation for the quasi- 
probability Q. Using the definition of the coherent 
states we find (see appendix |A| that equation to involve 
only first and second derivative terms respectively inter- 
pretable as classical drift and quantum diffusion, 



Q = (£ + C diS )Q 
C = id a u>a 



z + z* 



+ id z [-ujqz H — 7==(1 - z 2 ){a + a*)) + c.c. , £ di 



/2J 



(9) 

=d a d z (l - z 2 ) + c.c. , 



and may thus speak of a Fokker-Planck equation. We 
add in passing that we disregard any damping, restrict- 
ing ourselves, with respect to the condensate/cavity ex- 
periment of Ref . 6\ , to times smaller than the life times 
of both the cavity photons and the condensate. 

To discuss the semiclassical limit [13] j 3> 1, it is con- 
venient to switch to variables obeying canonical classical 
commutation relations. For the oscillator, we introduce 
'act ion- angle' variables (I, ip) through [H] 



a = V 3 Ic 



jle 



-lip 



(10) 



To parametrize the Bloch sphere of the spin we employ 
(cos #,</>), cf. ([3]). Both pairs are canonical, with Poisson 
brackets {/, tp}- = {cos 9, (/)} = 1. A few details pertaining 
to the change of variables are given in Appendix [B] 

Expressed in terms of these variables, the drift oper- 
ator C assumes the form of a Liouvillian, C = —{h, }, 
with the effective Hamiltonian function 



(11) 



h = ujq cos 9 + ljI + gy 81 cos ip sin 9 cos 



obtained from the Hamilton operator ([I]) by substitut- 
ing J — > jl, a — » \[JI exp(i^)), and dividing out hj = L. 
The classical approximation dtQ — CQ = — {h, Q} to 
the evolution equation ^ describes a drift of the quasi- 
probability Q(I, ip, cos 9, 4>) along the classical trajecto- 
ries of the Hamiltonian flow. The latter are determined 



by the Hamiltonian equations of motion 

/= — v / &g\/lsin'ipsm9cos(l) (12) 

ip — —uj — V2g —j= cos ip sin 9 cos <f> 

<p = uq — V8gvl cos ip cot 9 cos ef> 

cos9 = V8gVl costp sin9 smcp . 

As befits classical Hamiltonian equations the quantum 
number j does not show up here. 

When the quantum diffusion operator £diff is written 
in terms of the above canonical pairs of variables, it ac- 
quires a pre-factor \ which is very small in the semiclassi- 
cal limit. That semiclassical smallness notwithstanding, 
quantum diffusion has an important smoothing effect on 
the quasi-probabilty Q, as we shall see presently. But 
first, we devote a thorough discussion to the drift. 



III. CLASSICAL DYNAMICS 



In this section, we will analyse the phase space flow 



according to the classical Hamiltonian equations (12 1, in 
regimes of integrable, chaotic, and mixed dynamics. A 
discussion of chaos in the system has been reported in a 
seminal paper by Brandes and Emary [5]. However, the 
Holstein-Primakoff bosons employed to represent the spin 
variables in that reference tend to obscure the large-scale 
phase-space structure of the problem, and notably the 
semiclassical limits j 00 and h — > at hj = L = const. 
For an insightful discussion of chaos in the equations of 
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FIG. 1. Chaotic trajectories of the classical Dicke model at 
g — 3g c and h ~ 150cjo. The red and the blue trajectory 
differ in a slight mismatch of the initial conditions. 



motion (12 1 (applied to an opto-mechanical setting) we 
also refer to Ref. [TS]. The primary objective of our clas- 
sical analysis is to set the stage for the discussion of the 
quantum ramifications of chaos. 



A. Qualitative picture 



The equations of motion ( 12 ) do not involve the quan- 
tum parameters h and j = L/h. Remarkably, the scaled 
Hamiltonian h is also independent of the classical angu- 
lar momentum L — jh, due to the particular scaling ( 10 ) 
of the oscillator variables. Put differently, the magni- 
tude of the classical angular momentum does not affect 
the dynamics and can be accommodated in a rescaling 
of variables. The dynamics then depends on the dimen- 
sionless measures for frequency, coupling, and energy 
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h/uj , 



(13) 



with g c = ^loujq/2 the critical coupling at which the 
Dicke model undergoes its transition to a superradiant 
phase. Without much loss of generality, we will assume 
comparable frequencies uj/uj = 0(1) throughout [16] . To 
obtain some intuition of the dynamics, let us express the 
Hamiltonian in terms of the unit-length angular momen- 



h 

U 



= l z + vl + 7 v 2ul cos ipl x , 



(14) 



and consider high energy e 3> 1 and strong coupling 
7 > 1. The time variation of the oscillator phase is given 
by ip = —uj — g-^/2/I cos ipl x — uj + For suf- 

ficiently large I, and times larger than the oscillation 
period ~ uj^ 1 of ip, the second term becomes negligi- 
ble, i.e. we may approximate cos'0 — cos(u;t + ip(0))- 
The nearly harmonic oscillation of cos^ is exemplified 
by the trajectory shown in Fig. [I] (bottom left inset). It 
implies that the model behaves, at the high energies un- 
der consideration, much like a system of three dynamical 
variables (7, 6, <p) subject to external harmonic driving 
~ cosi/' at a frequency uj. The spin dynamics is governed 
by fast precession of the angular momentum 1 around the 
instantaneous rotation axis fl = uj^(^f\j2vl cos ip, 0, 1) T . 
For 'typical' values of cos?/>, the precession frequency 

^> uj ex- 
ceeds the 'driving' frequency uj by far. This is visi- 
ble in the fast spinning of the variable 1 around the 
unit sphere shown in the top section of Fig. [T] Second, 
the precession axis is typically oriented in x-direction, 
O ~ ujQ r )\/2vl cos tpe x + 0(1/ V I), and during these pe- 
riods the angular momentum component l x is approxi- 
mately conserved. This latter fact has important conse- 
quences for the variation of our primary variable of in- 
terest, I. Over time intervals of nearly conserved l x , the 
equation of motion / = — ujqj\/2iJI sin ipl x can be triv- 
ially integrated to obtain the characteristic arcs visible 
in the bottom right panel of Fig. [I] For any particular 
energy set by the initial condition, the action variable 
varies between an upper and a lower bound (indicated 
by horizontal lines) calculated in Appendix [Cj 

At times t ~ (mr + 7t/2)/cj (the ~ indicates a jitter 
of the order 1/y/l) the regular pattern outlined above 
gets interrupted, when the phase cos(wt) becomes small 
enough for Q x ~ cos(wt) and £l z to be comparable. Dur- 
ing these short time spans, the angular momentum pro- 
cesses around a vector ft no longer aligned in x-direction 
to a new orientation (cf. the isolated arcs visible in the 
top part of Fig. [I]) Specifically, the ir-component l x 
changes to a new and essentially un-predictable value. 
After the time window of small cosij) has been left, l x 
is approximately conserved again and the near regular 
change of I re-commences, at a changed rate ~ l x . 

Summarizing, the system behaves as if influenced by 
a 'random number generator': at regular time steps 
t (n7r + 7r/2)/tj, a new value of l x is dialed up, and 
that value sets the rate at which I(t) changes during 
the consecutive time interval of duration tt/2uj. Chaos 
manifests itself in this process's sensitivity to initial con- 
ditions. In the Fig. [TJthis is exemplified in terms of two 
trajectories of slightly different initial value of the coor- 
dinate I. By comparison, Fig. [2] visualizes the profile of 
trajectories in the integrable realm. Notice the strongly 
reduced, and effectively periodic fluctuations of the ac- 
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FIG. 2. Chaotic trajectories of the classical Dicke model at 
g — 0.2g c and h ~ 15CM;)- The red and the blue trajectory 
differ in a slight mismatch of the initial conditions. JfrNote 
the scale on which / varies, tiny compared to the chaotic case 
of Fig. 

tion variable, and the lack of divergence of trajectories of 
different initial conditions. 



B. Crossover to chaos 

Even though the literature on the Dicke model is vast 
no conclusive treatment is available of the emergence of 
chaos as energy and coupling strength are varied. Filling 
that gap appears all the more desirable as much previous 
work is focussed on low energies where the flow cannot 
explore all of the spin sphere. The relative status of the 
superradiant phase transition at g — g c and the crossover 
from regular to chaotic behavior could thus not be reli- 
ably ascertained. 

Global chaos is prevalent at large energies e > 1 and 
strong coupling 7 > 1. In general, the dynamics is mixed, 
or, in limiting cases, integrable. To map out the regimes 
of different dynamical behavior, we separately consider 
the model at weak and high excitation energies. 

Low energy dynamics . — At weak coupling 7 < 1, 
the Dicke Hamiltonian possesses a stationary point of 
lowest energy eo = — 1 a t I = 0, cos 6 = — 1. In the 
energetic vicinity of this point, the dynamics is inte- 
grable. Here 'vicinity' means excitation energies Ae <C 1, 
where eujQ = 2ujq defines the maximum energy that can 



be accommodated by the spin. Integrability is visible 
in the Poincare sections shown in the first few panels 
of Fig. [3] Signatures of mixed dynamics become vis- 
ible upon approaching the critical value 7 = 1. At 
7 > 1, the ground state configuration shifts to an en- 
ergy e = —{l 2 + 7~ 2 )/2, which is now attained at 
two degenerate points = ^("f 2 — 7 _2 ),0/7r) and 

(cos = (— 7~ 2 ,7r/0); here the non-vanishing value 
of the action coordinate corresponds to a macroscopic 
photon number (a^a) = jl. In the immediate vicinity 
of these points, the dynamics remains integrable, and 
for moderate excitation it is mixed. The subsequent 
crossover to chaotic dynamics turns out to be rather 
swift; already at excitation energies Ae ~ 0.2|e |, the 
energy shell is filled by chaotic trajectories. 

High energy dynamics . — At large energies, t > 1, the 
Bloch sphere gets fully covered by trajectories. Already 
at coupling strengths 7 < 1, trajectories become chaotic. 
The last tori get lost in the immediate vicinity of the 
critical value 7 = 1, cf. Fig. |4j 

For fully developped chaos, typical trajectories explore 
all of the energy shell eujQ = h. In view of our later discus- 
sion of photon number fluctuations, we need to explore 
the confines of these shells, especially with regard to the 
coordinate I. Lower and upper bounds of the action 
variable can be obtained as a result of a straightforward 
calculation detailed in Appendix [C] For large excitation 
energies, the accessible window of /-values asymptotes to 

with corrections of 0(1/ e). We are, thus, facing a window 
of width oc y/e and center oc e. In Section [V] we will 
consider the ramifications of the ergodic filling of these 
windows in the quantum dynamics of the system. 

C. Evolution of Q under chaotic drift 

Focusing on global chaos from now on, we would like to 
clarify how the Q-function would evolve if quantum dif- 
fusion were neglected. The assumed largeness of j indeed 
suggests to try out that most radical implementation of 
the semiclassical limit. 

For an initial coherent state, Q is isotropically 'sup- 
ported' by a single Planck cell, i.e., by a tiny fraction of 
the energy shell located at (o>q, Zo\H \oto, zq) = hui^e cx j 
with width hw Ae = [(ao,z \(H — hu> e) 2 \aQ, zq}] 1 ^ 2 cx 
y/j. By Liouville's theorem, that tiny fraction will not 
change in time. The chaotic evolution will interminably 
squeeze the originally 'circular support' of Q in the stable 
direction of the Hamiltonian flow and stretch it along the 
unstable direction. The ever narrower and longer 'sup- 
porting stripe' must soon begin to fold since the energy 
shell is compact. After a time of the order of the Ehren- 
fest time the stripe will have fully explored the energy 
shell. As the squeezing/stretching/folding of the support 
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FIG. 3. Poincare sections generated by monitoring the pair (l x , l y ) at fixed values of the phase ij), and at energies Ae = 0.2|eo| 
above the ground state; only southern hemisphere shown since northern one remains empty. For each parameter value, 7, nine 
trajectories of different on-shell initial conditions are sampled, a) 7 = 0.2, b) 7 = 0.7, c) 7 = 0.9, d) 7 = 1.01, e) 7 = 1.5. What 
looks like crossings in (b, c) is a reflection of the non-uniqueness of I(e, ip, cos 0, <j>), see App. B. Noteworthy is the predominance 
of chaos in (e), for an energy as small as 20% of the maximum energy capacity of the spin. 




FIG. 4. Poincare sections as in Fig. |3j now sampled at large energies Ae = 30|eo|. Values of 7 as in previous figure. The 
projection of the northern hemisphere, now fully covered by trajectories, looks qualitatively similar. 



of Q continues, an ever finer, and eventually 'singular' 
structure arises, where Q alternates infinitly rapidly be- 
tween high and near vanishing values transverse to the 
supporting stripe. Inasmuch as no region within the en- 
ergy shell appears favored, a constant mean density will 
arise. If one were to look at the 'fissured landscape' 
formed by Q with finite resolution one would, from a 
certain time on, just observe 'flatness' at the mean value 
of Q mentioned. In other words, one would see the micro- 
canonical distribution: Q constant within and zero out- 
side the energy shell. Expectation values of observables 
like low-order powers of the photon number ((a' a) m ) will 
not register the ruggedness of Q but just 'pick up' the 
microcanonical shape. Somewhat cavalierly said, Q ef- 
fectively equilibrates to the microcanonical distribution, 
within a time of the order of the Ehrcnfest time. 

The foregoing scenario changes little if we imagine the 
initial coherent state replaced by a squeezed minimum- 
uncertainty state. Initial states with larger uncertainties 
bring but two changes: (i) equilibration will happen even 
faster, the time scale shrinking logarithmically with the 
initial width, and (ii) the landscape underlying the effec- 
tively microcanonical Q can be smoother. 



IV. QUANTUM DIFFUSION 

Still confining ourselves to global chaos we now proceed 
to studying how quantum diffusion changes the effective 



cqulibration just found for the classical drift. We shall 
find a smoothing effect of quantum fluctuations which 
becomes effective, roughly, at phase space length scales 
~ y/h. To the best of our knowledge, our analysis of 
the Dicke system represents the first case study where 
the interplay of quantum fluctuations and nonlinear dy- 
namics in the long time behavior of a chaotic quantum 
system is resolved in concrete terms. A glance at the 
quantum diffusion operator (J9J) reveals that £diff couples 
oscillator variables to spin variables but does not include 
second-order derivatives wrt only oscillator variables nor 
wrt only spin variables. That structure is of course pre- 
served when the canonical pairs (I, ip) and (cos#, <j>) are 
employed, as we imagine done here. A real symmetric 
4x4 diffusion matrix then arises which has vanishing 'di- 
agonal' 2x2 blocks and mutually Hermitian conjugate 

'off-diagonal' 2x2 blocks, D = f° t fj . We shall not 

need the explicit dependence of the off-diagonal blocks 
d, d) on the variables /, ip cos 9, </> here but would like to 
emphasize the smallness d oc i . 

The 'chiral' block structure of D entails a secular equa- 
tion for the eigenvalues of the form A 4 — A 2 tr dd' + 
det d<V = 0. The four eigenvalues of D thus come in two 
plus/minus pairs ±Di,±D 2 where Df, D\ are the eigen- 
values of the non-negative 2x2 matrix dd) . Each of these 
pairs is associated with an eigenvector pair defining a 
contractive (— ) resp. expansive (+) direction. In the ex- 
pansive directions we confront normal diffusion while for 
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the contracting directions we may speak of anti-diffusion. but such cannot be expected any earlier than a Heisen- 

berg time tn oc j, possibly even a Poincare time (<x e J ). 



A. Qualitative discussion of equilibration 

The quantum diffusive contraction (expansion) com- 
petes, given chaos, with the stretching (shrinking) inher- 
ent in the classical drift. Quantum antidiffusive shrink- 
ing will be overwhelmed by the exponential determinis- 
tic expansion in the classically unstable direction: The 
pertinent scales will keep growing and the corresponding 
structure will ever more ubiquitously explore the energy 
shell, much as if quantum diffusion were entirely absent. 
However, in the deterministically stable direction where 
exponential shrinking proceeds ever more slowly, quan- 
tum diffusion will not allow that shrinking to go below 
a quantum scale oc Therefore, the Q-function will 
be diffusively smoothed transverse to the unstable di- 
rection such that the deterministically favored fissured 
'landscape' never arises. Microcanonical flatness will be 
reached from any initial state, coherent, squeezed, or 
broader, on the Ehrenfest time scale or faster. 

The foregoing arguments do not rule out revival events, 



B. Quantitative discussion: Co-moving quantum 
fluctuations 

The picture just drawn faithfully reflects a systematic 
theory obtained by projecting the dynamics onto a co- 
moving Poincare section. To see that we pick a phase 
space point X and a deflection 8X. The latter may 
be expressed through increments of the canonical vari- 
ables, SX = (SI, Sif), 5cos9, 5(j)) T , or through components 
(s,w, e, t) t = £ along the four directions distinguished 
asymptotically by the classical flow: one stable, one un- 
stable, and two neutrals (transverse to the energy shell 
and along the flow). The two variants are linearly related, 

(Sp, Sx, S cos 9, 8cj)) T — C(s, u, e, t) t , (16) 

or in brief 5X = C£. The 4x4 matrix C is com- 
posed by the unit vectors e s ,e u ,e e ,e T along the sta- 
ble/unstable/neutral directions as C = (e s , e Ul e e , e r ); it 
depends on Xq and must be determined numerically. 




FIG. 5. (Color online) Red: the coefficient D ss , Eq. numerically computed along a trajectory of energy h ~ 250u;o 

and a — 1.1. Blue: D ss integrated against the Lyapunov kernel exp(— 2Ai) as in (19 1. In the (arbitrary) units of the Plot, 
A" 



150. Dashed: the 'clock-variable' cos(f/>) effectively measuring time in units of oj 



Now we turn to the Fokker-Planck equation for Q. 
Near Xo the components of SX can be employed as 
Q = (d 5X A(Xo + SX) + dsxAx^aiXo + 5X))Q 
with d(X) the drift 'vector' and D(X) the diffusion ma- 
trix. We then change variables according to ( fl~6| and 
dsx = C T d^. The drift can be linearized and has the 
first two components d s = As, d u — —Am with A the 
Lyapounov rate. On the other hand, in the diffusion 
matrix we may drop the deflection £. Then integrating 
out the variables e, r we get the bivariate density Q(s, u) 
which obeys the reduced Fokker-Planck equation 

Q = (\d s s-\d u u+d 2 s D ss + dlD uu +2d s d u D su )Q . (17) 

Here the 'reduced diffusion matrix' matrix (^p ss J , 

describing quantum diffusion in the stable/ unstable sub- 
space, is the upper left 2x2 block in C T DC. Intuitive 



expressions arise when the orthonormal eigenvectors 
and eigenvalues of D are used, like 

D ss =^((e 5 ,V) 2 ^- (18) 
ft 

The sign of that effective diffusion constant will be of 
upmost importance. We see that the relative orientation 
of the intervening vectors matters, as well as the presence 
of positive and negative eigenvalues D^. 

Next, we proceed from the local quantum fluctuations 
near Xo to 'co-moving' fluctuations, simply following the 
classical trajectory starting at X . Save for the replace- 
ment Xo — > X t we have the same Fokker-Planck equa- 
tion at all times t, except that the 2x2 diffusion matrix 
becomes time dependent along the classical trajectory. 
The Lyapounov rate, on the other hand, neither changes 
along the trajectory nor when the point Xq is varied to 



select other (infinite) trajectories [T7] . 

The variance var t (s) = (s 2 ) t — (s) 2 of the stable deflec- 
tion is readily obtained as 



var t (s) 



e~ 2A *var ( 



{s) + / dt'e' 2X ^2D ss {t')\ (19) 



it must be positive at all times due to the guaranteed 
existence and positivity of the Q-function. In fact, to 
make sure we don't fall victims to our love of poetry, we 
have numerically checked that the diagonal element D ss 
remains mostly positive for a large number of trajectories 
of varying energy and coupling parameter. An exemplary 
plot of D ss evaluated along one of these trajectories is 
shown in Fig. [5] The plot exemplifies how D ss only rarely 
turns negative. The temporal convolution with e _2At in 
( 19 1 always entails a positive variance; see Fig. [5j 

Most remarkably, the local quantum fluctuations, man- 
ifest in the directions and strengths of diffusion and an- 
tidiffusion, 'sniff out' the asymptotically stable direction 
e s of the classical flow, tuning temselves to make for a 
positive variance varj(s) and thus a lower quantum bound 
~ 4^ for the scales accessible to the stable coordinate s. 
Co-motion is crucial since it lets the quantum fluctua- 
tions probe a time span at least of the order of the Lya- 
pounov time, over which stability properties of the classi- 
cal flow become manifest. Reassuringly, the linearization 
used to capture the co-moving fluctuations also remains 
reliable over that time span. 

In the unstable direction, on the other hand, the vari- 
ance var t (u) keeps growing indefinitely, as is similarly im- 
plied by the reduced Fokker-Planck equation (17). This 



is how the Q-function gets smoothed in the classically 
stable direction while forever extending its support in 
the unstable direction. The notorious singular structures 
of classical chaos are thus avoided, and equilibration to 
the microcanonical distribution takes place. 



V. GIANT FLUCTUATIONS DUE TO CHAOS 

Once more focusing on a range of energy and coupling 
strength with predominant classical chaos, we now turn 
to revealing giant fluctuations of the photon number in 
the stationary microcanical regime. 

If canonical pairs of variables are employed as argu- 
ments our equilibrated Q has the microcanonical form 



Q = ^S(h - ewo) 



(20) 



with ft the volume of the energy shell and e the (di- 
mensionless) energy set by the initial state. That micro- 
canonical equilibrium will of course be reached from any 
smooth initial state with an energy uncertainty similarly 
negligible as for a coherent state (where ^ ~ -^j). Writ- 
ing the delta function in the above distribution amounts 
to discarding corrections of relative order Station- 
ary means of powers of the photon number ((a^a) m ) for 



any integer m are now accessible as moments of the mi- 
crocanonical density ( 20 ) . In our semiclassical situation, 



these means are independent of the ordering of the 2m 
factors a,cv, to leading order in j. The normally or- 



dered form ((a 



measurable by absorbtion of the 



cavity output, is thus not different from the antinormally 
ordered form given by the moments of Q according to 
([6]) nor from the mean powers of the photon number, 
([a^a)™), such that we have 



(( a t ) 



dl /#<fcos0#r 



5{h — sujq) 



(21) 



The fourfold integral is most easily done in the case of 
high energies, e > 1, where the whole Bloch sphere is 
accessible (cf. Figs. 3j4 ) . As detailed in App. |Dj mean 
and variance of the photon number then come out as 



{at a) = (3 



((a f a) 2 ) - (a^a) 2 



UJ 

in 

3 la; 



euj 



3a; 
4ew g 2 



(22) 



136g 4 
15a; 2 



With the variance of order j 2 we indeed confront macro- 
scopic fluctuations. Notice that to leading order, and up 
to numerical factors, the estimate ( 15 ) is confirmed. 



We note without presenting calculations that the vari- 
ance of order j 2 persists down to smaller energies, pro- 
vided the energy shell is dominated by chaos. As visible 
in Fig. 2, e 0.2 suffices, together with 7 w 1.5. 



VI. BEYOND THE DICKE MODEL 

Evolution equations with derivatives terminating at 
second order are not restricted to the Dicke model. 
Whenever chaos is generated by a Hamiltonian of the 
form of a second-order polynomial in the pertinent ob- 
servables and a coherent-state-based Q-function can be 
used, we expect a Fokker-Planck equation for Q. Exam- 
ples are (i) SU(3) dynamics like the Lipkin model |18j . 
(ii) genuine many-body systems among which the Bosc- 
Hubbard model [HI is of much current interest (here 
the Hamiltonian is quartic in annihilation and creation 
operators, but due to the absence of antiresonant terms 
only first and second derivatives appear in the evolution 
equation for Q), and (iii) kicked systems like the top [T7] 
whose near classical quantum behavior has recently been 
observed experimentally |H] (see Appendix [E| . 

Even though numerous dynamical systems have 
Fokker-Planck equations representing their unitary quan- 
tum evolution, this behavior is by no means generic. 
In general, the Q-function evolves with derivatives be- 
yond the second order. For non-polynomial Hamiltoni- 
ans even infinite-order derivatives appear. The question 
then arises whether, given classical chaos, other equilibra- 
tion mechanisms reign or whether derivatives with orders 
n > 2 give but unimportant corrections to the quantum 
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diffusion carried by n = 2. Audacious as generai state- 
ments may be we dare pointing to a power counting ar- 
gument which suggests prevalence of the mechanism dis- 
cussed in this paper. When canonical pairs of variables 
are used, the generator of the time evolution of Q has the 
orders of derivatives and of Planck's constant interrelated 
as Yl n =i 2 ^ n l ^xfn[X) with X and the coefficients 
f n {X) independent of H. Herein n = 1 captures the clas- 
sical Hamiltonian drift while n — 2 accounts for quantum 
diffusion and brings about the minimal scale yh for the 
stable coordinates s. We may then set X — > X t in the 
coefficients /„ for n > 2, integrate out all but the sta- 
ble variables, and refer the stable variables to the said 
quantum scale as s = \Zhs. A reduced generator appears 
as dsXs + Y^n^ d™h( n ~ 2 ' > / 2 f n (X t ) and indeed suggests 
that quantum effects are dominated by the second-order 
derivative terms. 

Finally, inasmuch as a homogeneously filled energy 
shell has macroscopic extent in at least one phase space 
'direction', observables exploring that direction will dis- 
play macroscopic stationary fluctuations. 



VII. SUMMARY AND DISCUSSION 

The smoothing effect of diffusion on chaotic dynamics 
has been noted before, e.g., within the context of quan- 
tum billiards (cf. Ref. 22J). However, diffusive contri- 
butions to classical evolution were there added by hand. 
Our present analysis exemplifies how unitary quantum 
evolution itself brings about diffusion. By projecting the 
Dicke model dynamics onto a co-moving Poincare surface 
of section we could check explicitly that quantum diffu- 
sion sets a limiting scale to the variance of the stable 
coordinate such that Q equilibrates to a smooth density 
of the microcanonical form. 

It is well to realize that we are facing a privilege of the 
Q-function which other popular quasi-probability densi- 
ties like the Wigner function or the Glauber-Sudarshan 
P-function (weight in a diagonal mixture of coherent 
states) do not enjoy. The Wigner function W, for in- 
stance, is known to develop positive/negative substruc- 
tures within Planck cells under conditions of classical 
chaos 23 . Such substructures forbid pointwise conver- 
gence of W to a classical probability density as h — » 0; 
they are washed out by the average over, roughly, a 
Planck cell which leads from W to Q. The situation 
is even more precarious for the P-function from which 
W arises by smoothing over, roughly, a Planck cell. Not 
only is P prone to going negative but even to loosing 
existence as an ordinary function under dynamics with 
classical chaos. For instance, a coherent initial state will 
get its support distorted to that of a 'Schrodinger cat 
state' (in the classicaly unstable direction) [23], which 
latter is known to have a non-positive and even singular 
P 24J . It is in fact easy to check that for the Dicke model 
the diffusion terms Ajiff for P and Q differ only in sign 
(see Appendix |A| ; therefore, the variance var t (s) which 



remains positive at all times for Q must sooner or later 
go negative for P. 

In response to the recent experimental observation of 
the supcrradiant phase transition in the Dicke model, 
we have investigated the prospects of detecting the con- 
comitant transition from regular dynamics at the low- 
est of energies to prevalence of chaos at higher excita- 
tions. As a most interesting witness of that transition 
we have identified stationary fluctuations of the number 
of oscillator quanta (photons). While small for regular 
dynamics, these fluctuations rise to macroscopic magni- 
tude as chaos proceeds towards fully covering the energy 
shell. Perhaps fortunately for attempts at detection, the 
large-fluctuation regime signalling fully chaotic behavior 
is found already for moderate degrees of excitation, pro- 
vided the coupling is chosen above the critical value for 
the superradiant phase transition. The giant fluctuations 
are predicted to arise independent of the initial state, af- 
ter a time of the order of the Ehrcnfest time. 

We have also argued that both our equilibration mech- 
anism and large fluctuations of suitable observables are 
at work in other observable systems of current interest. 
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Forschungsgemeinschaft . 



Appendix A: Derivation of the quantum evolution 
equation 

We start the derivation of Eq. ([9| with the identity 
2j + 1 

Q( a > z ) = ^2(1 + U|2)2 \ Q ' z \ d tp\a, z) 

= ~ L 2 a +^zp) 2 tr ^' |a ' z)(a ' 2:| ] - (A1) 

To process this expression, we need to compute the action 
of the Hilbert operators on coherent states. As a result 
of a straightforward calculation, one obtains 

J-\z){z\ = (d z + 2j I ^\ \z)(z\, 

J+\z){z\ = (-z*d z + 2j^-^ \z){z\, 

l\z){z\ = {-zd z l*M (A2) 

and 

a|a)(a| = a|a)(a|, 

a?\a){a\ = (d a + a*)\a){a\. (A3) 
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Substituting Eq. ([I]) into (Al) and using the relations 
above we obtain our Fokker-Planck equation 

Likewise, one checks that the Glauber-Sudarshan P- 
function, defined as the weight in the diagonal mixture 
of coherent states 



p= I dadzP(a, z)\a, z) (a, z\ , 



(A4) 



obeys a Fokker-Planck equation whose generator differs 
from the one for the Q- function (cf ([9])) in only two de- 
tails: (i) The factor (j + 1) in the drift term in ^ is 
replaced by j. (ii) Much more importantly, the diffusion 
term acquires an overall minus sign. 



D = 



and the diffusion matrix 
d\ 

Die Die/, 



(B5) 



l 

l-c 2 



A = — cos ip sin <j) — c sin ip cos <f> , 
B = c cos ip cos (j> — sin ip sin <p . 



Appendix C: Bounds on photon number in energy 
shell 



Appendix B: Transformation to canonical variables 

For completeness, we here provide a few technical de- 
tails relating to the change of variables (a,a* , z, z*) — > 
(I, ip,c= cos 9, </>). The defining relations 



z = tan-e i *= J— c" 
2 V 1 + c 



yield the derivatives 



jle 



lip 



(Bl) 



d c ( i + c) + id^ \j\^ c 



d a = 




(B2) 



-ii/j 



and their complex conjugates. We must realize that 
the complex stereographic projection variables (z, z*) are 
not a canonical pair. Therefore, a Jacobian arises in 
Q(c,^,I,ip) = ■^^Q(z(c,<l>),z*(c ) <f>),a(I > 'ip),a*(I,i))) 

such that to get the generator for Q(c,<f>,I,if)) we must 
replace as d c -> (l + c)~ 2 9 c (l + c) 2 = <9 c + 2(l + cr 1 and 



9 c (l + c)Vl-c 2 +i5 v 



1 + c 



l-c 



(B3) 



Straightforward calculation then gives the generator in 
search as dxidi + dx^XjOij with the drift vector 



di 
d c 



y/8g(l + — ^ VTsin i\)\J 1 — c 2 cos 
cj + \/2#(l + cosVVl - c 2 



(B4) 



-a>o + V$>gVl cos ip 



-V&gVl cos ipyjl — c 2 sin^, 



Solving Eq. (14) for /, one readily finds 
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±W (ccsin( 



(e - cos 6>) - x sin 9 , (CI) 



where we defined x = cos ip cos 0. Interestingly, as a func- 
tion of e, cos 0, (/>, ?/;, the quantity / may take on two dif- 
ferent values, and that complication arises when both x 
and e — cos 9 are negative. No solution for / exists when 
non-negative x meets with negative e — cos 9. The sim- 
plest situation is e — cos 9 > 0: then only the positive 
square root in (CI) is possible and / is unique. We illus- 



trate the search for the bounds with just a few cursory 
remarks on the latter case. 

There is no extremum in the calculus sense. So the 
smallest and largest values of / must occur on the bound- 
aries cos 9 = ±1 or/and x = ±1. The poles of the Bloch 
sphere provide the "trivial" bounds /mm = i/ _1 (e ± 1). 

To check the possibility of tighter bounds we first try 
x = 1. From d\fl / d cos 9 = we get 



cos 9 



7 2 ± ^7 4 + 1 + 2 7 2 
1 + 2 7 2 e 



7 



'2v 



tan0. (C2) 



Due to \fl > and the global positivity of sin 0, we 
must require cos 9 > and therefore only the upper sign 
qualifies. A non-trivial lower bound / m i n is thus obtained. 
Similarly, the case x = — 1 yields a non-trivial upper 
bound v 7 /max — — tan 9 with tan 9 according to the 



lower sign in (C2). The leading terms of the e-expansion 



of these bounds are the ones given in ( 15 ). 



Appendix D: Microcanonical averages 

We briefly sketch the calculation of the microcanonical 
moments (21), for simplicity confining ourselves to high 



energies (e > 1). Doing the /-integral in (21) we have 



1 f I m 
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with the peak intensity / = I(ip, <j>, 9) determined by ( CI ) 



and x = cos t/> cos 4>; only the positive sign in ( CI I is pos- 



sible at high energies. The derivative of the Hamiltonian 
1 2 I 



\dh/dl\ ujoV^I ^J( X sm9) 2 + ^{e-cos9) 



allows to rewrite the moments as 

^2 ^ m i 



M„ 



rr 

1v 



dip dcos 9 d(j) 



( 



JJx sin 6) 2 + |r (e - cos 9) 



x sin ( 



2 m +1 



JJx sin 9) 2 + ^ (e - cos 9) 



The power in the foregoing numerator can be binomially 
expanded. By symmetry only even powers of x sin 9 con- 
tribute and therefore only even powers of the square root 
remain. We quickly find 



Mi = 

M 2 = 



(< 2 + + i-,' + J) 



which immediately implies ( 22 1 . 



Appendix E: Kicked top 

We would like to corroborate our expectation for the 
kicked top [181 125] , a periodically kicked large spin with 

conserved length, J 2 = j(j : + 1) 3> 1. Classical equilibra- 
tion for a cloud of points on the Bloch sphere has long 
been known from numerical studies. Quantum equilibra- 
tion and the ensuing large stationary fluctuations of the 
(orientation of the) angular momntum should be observ- 
able in a variant of the experiment of Ref. [21] . 
The simplest chaotic top has the Floquet operator 



F 



it involves a rotation about the J^-axis by the angle p and 
a subsequent 'torsion' about the J 2 -axis. Torsion means 
a state dependent rotation by the angle ^J+l which has 
opposite signs in the northern and southern hemispheres. 
The precession angle p and the torsion constant r are 
assumed independent of j. Chaos predominates if r 3> 1. 

The stroboscopic time evolution of the density oper- 
ator is given by p n = F n poF' n with the dimension- 
less integer 'time' n. Employing the Q-function Q(z) = 
■k{i+zI-) 2 we §° f° r the propagator for its single- 

step evolution Q n +\{z) = FQ n {z). Like the Floquet op- 
erator F, the Husimi propagator J- = J- T J- p is a product 
of two factors, one each for precession and torsion. 

For the precession we get 

T p = cxp {^ipd z (l - z 1 ) + c.c.} . 
The generator in the foregoing exponent involves only 
drift (first order derivative terms) but no diffusion; it 
is the generator for rotation about the J x -axis already 
encountered for the Dicke model in classical Hamil- 
tonian in character. The torsion propagator reads 



T T = cxp< —it 



1 — zz* 
1 + zz* 



2j + l 



d 2 z 2 + c.c 



with an exponent involving drift and diffusion. A i- 
correction in the drift has not been written out. The 
displayed drift is again classical Hamiltonian, as becomes 
visible once the real canonical pair cos 9, <f> of variables is 
introduced. The quantum diffusion with an explicit fac- 
tor of order \ involves a real symmetric 2x2 diffusion 
matrix with vanishing diagonal elements (i. e., of chiral 
structure). 

All arguments for equilibration for the Dicke model ap- 
ply again. Most importantly, the positive eigenvalue of 
the diffusion matrix sets a smallest scale for the motion 
along the stable direction of the classical drift. On that 
latter scale Q becomes smooth transverse to the classi- 
cally unstable direction. Effective stationarity will reign 
no later than about an Ehrenfest time, with Q constant 
over the Bloch sphere. Equipartition of Q will result in 
large stationary fluctuations of the angular momentum. 
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